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Abstract 

We present a next-to-leading-order calculation of the production of a Z boson with 
two jets, one or more of which contains a heavy quark {Q = c,b). We show that 
the cross section with only one heavy-quark jet is larger than that with two heavy- 
quark jets at both the Fermilab Tevatron and the CERN LHC. These processes are 
the dominant irreducible backgrounds to a Higgs boson produced in association with 
a Z boson, followed hy h ^ bb. Our calculation makes use of a heavy-quark distribu- 
tion function, which resums collinear logarithms and makes the next-to-leading-order 
calculation tractable. 



1 Introduction 



The discovery of new physics at hadron coUiders often rehes on a detailed understanding of 
standard-model background processes. Prominent among these is the production of weak 
bosons {W, Z) in association with jets, one or more of which contains a heavy quark {Q = 
c, b). The prime example is the discovery of the top quark at the Fermilab Tevatron, which 
required a thorough understanding of the W^+jets background, with one or more heavy-quark 
jets [HEIE]. The discovery of single-top-quark production via the weak interaction requires 
an even more sophisticated understanding of this background [H [7] . 

In this paper we present a next-to-leading-order (NLO) calculation of the production 
of a Z boson in association with two or more jets, one or more of which contains a heavy 
quark. This is the dominant irreducible background to the production of a Higgs boson in 
association with a Z boson, followed hj h ^ bb [HI E]. This signal is currently being sought 
at the Tevatron with either one or both b jets tagged [8] . The calculation we present provides 
the irreducible background at NLO for both cases. 

The production of a Z boson plus two jets with at least one heavy-quark jet is also a 
background to the production of a Higgs boson in association with one or more b jets, which 
is a discovery mode for a supersymmetric Higgs boson at large values of tan j3 [H JTUl ttH 

[I21[I3l[Il[iaESltt3ttSl[ia[2ni[ni[221[23l[MlESlESl[2^ EgHrhis 

background process is also a benchmark for this Higgs discovery channel. The search for a 
supersymmetric Higgs boson via this mode is underway at the Tevatron [33l [311 ES] , and will 
be vigorously pursued at the CERN Large Hadron Collider (LHC) [361 EZ] ■ 

When one considers the production of a heavy quark at a hadron collider, ones first 
thought is usually of a virtual gluon splitting into a final-state QQ pair, as shown in Fig.[l](a), 
or via gluon fusion, as shown in Fig. [T](b). However, initial gluons splitting into a QQ pair 
is just as important a source at the Tevatron, and even more important at the LHC, when 
only one heavy quark is observed at high transverse momentum (pt) ■ In that situation, it is 
advantageous to think of the initial gluon as splitting into a collinear QQ pair, with one heavy 
quark remaining at low pt while the other heavy quark participates in the hard scattering and 
emerges at high pt- The heavy quark that participates in the hard scattering can be treated 
as part of the proton sea, with a parton distribution function that is calculated perturbatively 
from the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [3, [6]. 
The production of a Z boson with two jets, one of which contains a heavy quark, then 
proceeds as shown in Fig. [21 

The reasons for using a heavy-quark distribution function are twofold. First, it resums 
collinear logarithms of the form In Q/itiq to all orders, where Q is the scale of the hard scat- 
tering. Second, it simplifies the leading-order process, which makes a higher-order calculation 
tractable. 

This paper completes the NLO calculation of Z+jets, with one or more heavy quarks, 
up to two jets [IQ]. The production of a Z plus one jet, with one or more heavy quarks, 
was presented in Ref. [H], and agrees well with data from the Tevatron [12]. The inclusive 
production of a Z with one or more heavy quarks was presented in Ref. [33] • 

^The minimal supersymmetric standard model requires two Higgs doublets; the ratio of their vacuum 
expectation values is tan f3 = V2 /vi . 
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Figure 2: Diagrams contributing to the associated production of a Z boson and two high-py 
jets, one of which contains a heavy quark {Q = c,b). 

2 ZQj at NLO 

The leading-order (LO) processes for Z plus two jets, one or more of which contains a heavy 
quark, are qq{gg) — > ZQQ (Fig. [1]) and Qq{g) — * ZQq{g) (Fig. [2]). The LO cross sections 
are given in Table [1], with the jets satisfying the conditions pt > 15 GeV, \t]\ < 2 (2.5 
at the LHC), and ARjj > 0.7. At the Tevatron, the leading-order cross sections for these 
two classes of processes are comparable; at the LHC, the latter dominates. The processes 
Qq{g) ~^ ZQq{g) relatively more important at the LHC than at the Tevatron because 
they are initiated by a heavy sea quark, whose distribution function rises at small values of x. 
Furthermore, Qg ZQg is larger than Qq ZQq at the LHC, while they are comparable 
at the Tevatron, since the former involves the gluon distribution function, which is large at 
small values of x. For the same reason, gg —>■ ZQQ is dominant at the LHC compared with 
qq —>■ ZQQ, while the latter is more important at the Tevatron. 

As is often the case, the distinction between the various processes is lost once one goes 
beyond leading order. Fig. [2](a) shows a Feynman diagram which ostensibly contributes 
to the NLO correction to qq ZQQ, while the Feynman diagram in Fig. Mjo) ostensibly 
contributes to the NLO correction to Qq ZQq, and the diagram in Fig. [2](c) to the NLO 
correction to gg ZQQ. However, all three diagrams contribute to the same amplitude, and 
therefore interfere. Thus one cannot uniquely identify them with any of the leading-order 
processes. 

The next-to-leading-order calculations in this paper were performed with the Monte-Carlo 
code MCFM |36]. The leading-order calculations were performed both with this code and 
with MadEvent [TTj. The NLO corrections are included in the code MCFM by implementing 
the virtual helicity amplitudes of Ref. [58]. These are given in the four- dimensional helicity 
scheme, which is used throughout the calculation. The real corrections are adapted from 
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Table 1: Leading-order cross sections (pb) for Z boson plus two jets, one or two of which 
contains a heavy quark, at the Tevatron (-^5 = 1.96 TeV pp) and the LHC (a/s = 14 TeV 
pp). A jet lies in the range p-r > 15 GeV and \!]\ < 2 (2.5 at the LHC), with ARjj > 0.7. 
No branching ratios or tagging efficiencies are included. The labels on the columns have 
the following meaning: ZQj = exactly two jets, one of which contains a heavy quark; 
ZQQ = exactly two jets, both of which contain a heavy quark. The CTEQ6L1 parton 
distribution functions are used with the factorization and renormalization scales chosen 
as /ii? = /ii? = Mz- Also given, in square brackets, is the LO cross section obtained without 
using a heavy-quark distribution function, from gq{g) ZQQq{g). 



Process 


a(pb) 


Tevatron 


LHC 


ZQj ZQQ 


ZQj ZQQ 


bq Zbq + bg —>■ Zbg 
qq — s> Zbb + gg ^ Zbb 


0.89+1.29=2.18 [1.91] 

1.89+0.58=2.47 


76+276=352 [191] 

13+96=109 


cq — s> Zcq + eg ^ Zcg 
qq — >• Zee + gg ^ Zee 


1.37+1.83=3.20 [3.26] 

1.89+0.45=2.34 


98+345=443 [271] 

12+75=87 
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(a) (b) (c) 

Figure 3: Diagrams contributing to the NLO correction to the associated production of a Z 
boson and two high-p^ jets, one or more of which contains a heavy quark {Q = c,b). 

Ref. jl^ and singularities are handled using the dipole subtraction method [1]. Since the 
matrix elements contain the decay of the gauge boson into massless particles, the code is 
general enough to provide results for final states from 7*, Z* (including interference). 

For this paper we specialize to the case of a real Z boson. 
The processes involved in the calculation are as follows: 
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We also include the processes QQ' ZQQ', QQ ZQQ , and QQ ZQ'Q at LO. These 
processes are already small at LO, so it is safe to neglect their NLO corrections. For charm 
final states, we only take Q = Q' = c; for bottom, Q,Q' = c,b. 

The results of our calculation are presented in Tables [2] and [31 The columns ZQj and 
ZQQ contain the leading-order cross sections (in parentheses) for the processes Qq{g) — * 
ZQq{g) and qq{gg) — > ZQQ, respectively, taken from Tabled! We require the jets to satisfy 
the conditions pt > 15 GeV, |?7| < 2 (2.5 at the LHC), and ARjj > 0.7. The heavy- 
quark mass is neglected here and throughout the calculation (except where noted). To 
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Table 2: Cross sections (pb) for Z-boson plus two (or more) jets, one or more of which 
contains a heavy quark, at the Tevatron {^/s = 1.96 TeV pp). A jet lies in the range pt > 15 
GeV and \7]\ < 2. Two final-state partons are merged into a single jet if ^Rjj < 0.7. No 
branching ratios or tagging efficiencies are included. Numbers in parenthesis are leading- 
order results. The labels on the columns have the following meaning: ZQj = exactly two 
jets, one of which contains a heavy quark; ZQQ = exactly two jets, both of which contain 
a heavy quark; Z{QQ)j = exactly two jets, one of which contains a heavy-quark pair; 
ZQQj = exactly three jets, two of which contain a heavy quark; ZQjj = exactly three jets, 
one of which contains a heavy quark. For the last set of processes, the labels mean: Zjj = 
exactly two jets, including heavy quarks; Zjjj = exactly three jets, including heavy quarks. 
For ZQj and ZQQ, both the leading-order (in parentheses) and next-to-leading-order cross 
sections are given. The CTEQ6M parton distribution functions are used throughout, except 
for the LO cross sections in parentheses, where CTEQ6L1 is used The factorization and 
renormalization scales are chosen as fip = fiR = Mz- 



Tevatron 


a (pb) 




ZQj ZQQ ZiQQ)j ZQQj ZQjj 


bottom 
charm 


(2.18) 5.23 (2.47) 3.07 0.634 0.672 0.326 
(3.20) 7.49 (2.34) 2.75 2.00 0.621 0.495 




Zjj Zjjj 


Z+jets 


(163) 182 22.9 



obtain the NLO cross sections for ZQj and ZQQ, which involve the radiation of additional 
partons {e.g.. Fig. [2]), two partons are combined into a single jet by adding their four- 
momenta if ^Rjj < 0.7. If the combined partons are both heavy quarks, then the process 
contributes to the column labeled Z{QQ)j. This is a Z -|- 2j event in which one jet contains 
two heavy quarks, which changes the tagging probability for that jet [51]. It is calculated 
with a finite heavy-quark mass {iric = 1.4 GeV, rrif, = 4.75 GeV) in order to regulate the 
logarithmic divergence present when the heavy quarks are coUinear. If all three partons are 
well separated, then the process contributes to either ZQQj or ZQjj^ 

We checked that the effect of the heavy-quark mass is negligible by comparing ZQQ at 
tree level with and without a finite quark mass. Similarly, we found that the heavy-quark 
mass is also negligible for ZQQj at tree level, as expected. 

The NLO correction to ZQj is quite sizable at the Tevatron, more than 100%. One 
of the reasons is that the process qq ZQQg (where one of the heavy quarks is outside 
the acceptance) makes a relatively large contribution. This is not really a NLO correction 

^The matrix elements for both Z{QQ)j and ZQQj are calculated with a finite heavy-quark mass, but 
massive phase space is used only in the former. We checked that the quark mass is numerically irrelevant to 
the latter process. 
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Table 3: Same as Table [2], except at the LHC {^/s = 14 TeV pp). A jet lies in the range 
Pt > 15 GeV and |?7| < 2.5. 



LHC 


a (pb) 




ZQj ZQQ Z{QQ)j ZQQj ZQjj 


bottom 
charm 


(352) 421 (109) 92.1 23.5 60.8 92.1 
(443) 623 (87) 75.2 58.6 49.3 123 




ZjJ Zjjj 


Z+jets 


(6090) 4840 1810 



to Qq{g) — > ZQq{g), but rather a new channel. This contribution is about 1.2 pb for 
both bottom and charm at the Tevatron. In contrast, this process makes a relatively small 
contribution to ZQj at the LHC, only 10 pb for bottom and 4 pb for charm. 

We also list in Tables [2] and [3] the LO (in parentheses) and NLO cross sections for Zjj 
[ini [S21 [S3], and the LO cross section for Zjjj. In these cross sections we have included 
the contribution from light partons as well as heavy quarks. Thus, for example, the fraction 
of Z + 2j events in which only one of the jets contains heavy quarks is given by [ZQj + 
ZiQQ)j]/Zjj. 

In Tables H] and [5] we list the results in an inclusive manner, and also include the NLO 
result for Z plus one heavy-quark jet ^1]. In addition, we give the NLO cross section for 
Zj [E], including both light and heavy partons. Thus ZQ + X is the cross section for a Z 
plus at least one heavy-quark jet; it is the sum of the columns labeled ZQ, ZQj, and ZQQ 
in Tables 1 and 2 of Ref. ^Ij, including all contributing processes. Z{QQ) contains one jet 
which contains a heavy-quark pair. It is calculated at LO, including the heavy-quark mass 
in order to regulate the collinear divergence present in qq Z{QQ). ZQj + X has at least 
two jets, one of which contains a heavy quark; it is the sum of the columns labeled ZQj 
and ZQjj in Tables |2] and [31 ZQQ + X has at least two jets, two of which contain a heavy 
quark; it is the sum of the columns labeled ZQQ and ZQQj in Tables [2] and [31 Finally, 
Z{QQ)j contains two jets, one of which contains a heavy-quark pair; it is calculated at LO 
with a finite quark mass. 

We estimate the uncertainties in the NLO inclusive cross sections by varying the renor- 
malization scale, the factorization scale, and the parton distribution functions independently. 
The renormalization scale is varied over = {0-5 — 2)Mz, with fip = Mz- Similarly, we vary 
fj,F = (0.5 — 2)Mz with /i/j = Mz- The uncertainties due to scale variation are significantly 
reduced at NLO in comparison with the LO calculation. Significant renormalization-scale 
dependence remains since the process is (9(a|). The parton distribution functions are varied 
over the 41 different sets contained in CTEQ6M [71 [IS]. A symmetric uncertainty on the 
cross sections can be obtained from these sets by using Eq. (26) of Ref. [53] . In addition, we 
have studied the deviation of the cross section in each direction separately, as in Ref. [55] . 



6 




50 100 150 200 

Pt(Z) [GeV] 

Figure 4: Transverse-momentum distribution of the Z boson in events with at least two jets, 
one of which contains a bottom quark, at the Tevatron (-y/i = 1.96 TeV pp). The jets satisfy 
the conditions -pr > 15 GeV, \r]\ < 2, and ARjj > 0.7. Both LO (dashed) and NLO (sohd) 
distributions are shown. 

However, we find that the difference between the positive and negative deviations is not 
significant in any of the cases and thus we report only symmetric uncertainties. 

As an example of the use of these tables, consider the inclusive cross section for Z+ 2 
jets with one heavy-quark tag. This cross section is obtained from 

O- = eqCTzQj+X + 2eQ(l - eQ)o-^QQ+x + ^QQ^Z{QQ)j (1) 

where eg is the tagging probability for a heavy-quark jet, and Cgg is the tagging probability 
for a jet containing a heavy-quark pair. 

We show in Figures [Hand [5] the transverse- momentum spectrum of the Z boson in events 
with at least two jets, one of which contains a bottom quark, at the Tevatron and the LHC 
Both the LO and NLO distributions are shown. The radiative corrections do not significantly 
change the shape of the distribution, as is often the case when the quantity plotted is not 
affected by the change in the kinematics due to additional radiation. 

3 Conclusions 

In this paper we present a NLO calculation of the production of a Z boson plus two jets, one 
or more of which contains a heavy quark. This greatly improves the accuracy with which 
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Table 4: Inclusive cross sections (pb) for Z boson plus one or two jets, one or two of which 
contains a heavy quark, at the Tevatron (y^ = 1.96 TeV pp). A jet lies in the range > 15 
GeV and |?7| < 2. Two final-state partons are merged into a single jet if Ai?jj < 0.7. No 
branching ratios or tagging efficiencies are included. Numbers in parenthesis are leading- 
order results. The labels on the columns have the following meaning: ZQ + X = at least 
one jet, at least one of which contains a heavy quark; Z{QQ) = one jet which contains 
a heavy-quark pair; ZQj + X = at least two jets, one of which contains a heavy quark; 
ZQQ + X = aX least two jets, two of which contain a heavy quark; Z{QQ)j = two jets, 
one of which contains a heavy-quark pair. The last row gives the inclusive cross section for 
jets containing both light and heavy partons. The CTEQ6M parton distribution functions 
are used throughout, except for the LO cross sections in parentheses, where CTEQ6L1 is 
used |7]. The factorization and renormalization scales are chosen as /i^ = /^k = Mz- The 
uncertainties are from the variation of the renormalization scale, the factorization scale, and 
the parton distribution functions, respectively. 



Tevatron 


a (pb) 




Z + 1 jet +X 


Z + 2 jets +X 


ZQ + X Z{QQ) 


ZQj + X ZQQ + X Z{QQ)j 


bottom 
charm 


(8.23) 18.1 2.1 
(11.3) 27.5 6.6 


(0 ^Q^ c; c;(5+1.2 +0.07 +0.5 ^9 /1q\ o 7/1+0.45 +0.12 +0.15 n fi'i+O ^e +0.06 +0.05 
I^Z.iyj 0.00_o.9 _o.05 -0.5 [^-^y) <J- '"±-0.45 -0.12 -0.15 U.0t)_o.l6 -0.06 -0.05 

/"? 91 S 9'i+l-'' +0.15 +0.8 /9 ocr\ /I 7+0.44 +0.0 +0.14 9 +0-22 +0.15 
yO.Zl) o.ZO_;^ 4 _o.26 -0.8 y^-'J'J) <J-^' -0.37 -0.87 -0.14 ^•Uo_o.53 -0.16 -0.15 


all jets 


(898) 1070 


(163) 205t\ltin 



Table 5: Same as Table HI except at the LHC {y/s = 14 TeV pp). A jet lies in the range 
Pt > 15 GeV and \t]\ < 2.5. 



LHC 


a (pb) 




Z + 1 jet +X 


Z + 2 jets +X 


ZQ + X Z{QQ) 


ZQ] + X ZQQ + X Z{QQ)3 


bottom 
charm 


(826) 1060 25 
(989) 1430 50 


CiKo-X r^^o+?>i+ii+25 / 1 i i n i p-q+20 +2 +9 9/1+10 +0.3 +1.2 
I^OOOJ 010_5g _35 _25 I 1 ij 100_20 -2 -9 ^^-6 -0.3-1.2 

(443) 746111° t% tt (90) 12511? li tl 59li tl H 


all jets 


(15300) 18400 


(6090) 6650l^^r™l^l° 
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Figure 5: Same as Fig. HJ but at the LHC {^/s = 14 TeV pp). The jets satisfy the conditions 
Pt > 15 GeV, |?7| < 2.5, and A% > 0.7. 

this fundamental background is known at the Tevatron and the LHC. We provide our results 
in both an exclusive (Tables [2] and [3]) and an inclusive manner (Tables H] and [5]). The NLO 
cross section is significantly greater than the LO cross section for the default scale choice 

IJ-F = fJ'R = Mz- 

Our calculation makes use of the heavy quarks present in the proton sea, which are 
perturbatively calculable. This makes the LO calculation simpler, thus allowing a higher- 
order calculation to be tractable. We showed that the processes Qq{g) — > ZQq[g) are a 
significant source of ZQj events at the Tevatron, and the dominant source at the LHC. 

Alternatively, one could eschew the heavy quarks in the proton sea, and regard the proton 
as containing only light quarks and gluons. In that approach, the LO processes for ZQj are 
gi{g) ~^ ZQQq[g)^ where one of the heavy quarks is outside the acceptance of the detector. 
We give the results of this calculation in square brackets in Table [T], using CTEQLl with 
fJ^F = fJ'R = Mz- The results are quite consistent with the LO results of the heavy-quark 
approach at the Tevatron, but at the LHC they lie somewhat below, though within a factor 
of two. In any case, the most accurate cross sections are the NLO results presented in this 
paper, based on the heavy-quark approach. 

An important application of these results is to the search for the Higgs boson via qq Zh, 
followed hj h ^ bb, where the signal is Z + 2j with one or two b tags [8j. We provide NLO 
results for the backgrounds with either one or two heavy-quark jets. We have shown that 
the majority of such events have only one heavy-quark jet. We see from Tables H] and [5] that 
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Zhj + X is nearly twice as big as Zhh + X at the Tevatron, and more than three times as 
large at the LHC Similarly, Zcj + X is about thrice Zee + X at the Tevatron, and about 
six times as large at the LHC. 
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We recently presented a next-to-leading-order (NLO) calculation of Z + 2 jets, with one 
or more heavy-quarks {Q = c,b), at hadron colliders [1] . This process is a background to 
numerous searches at hadron colliders, in particular the search for the Higgs boson: it is the 
principal irreducible background to qq — >■ Zh, followed by /i — > bb, and a significant reducible 
background to gg —>■ h ZZ — >■ when both jets contain heavy quarks that decay 

semi-leptonically. We neglected the heavy-quark mass throughout that calculation, except 
when two heavy quarks were produced from a virtual gluon, g* —>■ QQ, in which case it is 
mandatory to keep the quark mass nonzero when the heavy quarks are coUinear. This occurs 
when two heavy quarks are contained in the same jet. More generally, one cannot neglect 
the quark mass if the invariant mass of the quark pair coming from g* QQ is not large 
compared with the quark mass. 

In a subsequent NLO calculation of W + 2 jets, with one or more b quarks, we learned 
that this can also occur in another situation [2j, which we overlooked in Ref. [T|. If both 
heavy quarks are at large transverse momentum (pt) compared to their mass, then the mass 
may be safely neglected [3J. However, there are NLO processes in which one heavy quark 
is at high pt while another is at low pt (and is missed by the detector), yet their invariant 
mass is not large compared to their mass. If this heavy-quark pair comes from a virtual 
gluon sphtting to a heavy-quark pair, g* QQ, then one cannot neglect the quark mass. 
This can occur in the NLO processes qq ZQQg (Fig. [1]) and gq ZQQq (Fig. [2]), both 
of which contribute to ZQj when one heavy quark is missed. 



Figure 1: Diagram contributing to the NLO correction to the associated production of a Z 
boson and two high-pr jets, one or more of which contains a heavy quark {Q = c,b). 

To demonstrate this, we show in Fig.Othe invariant-mass distribution of the heavy-quark 
pair from qq Zbbg when it contributes to the final state Zbj, with the other heavy quark 
outside the acceptance of the detector (taken to be px > 15 GeV, \ri\ < 2, ARjj > 0.7), at 
the Fermilab Tevatron (v^ = 1.96 TeV pp). There are two curves, one with (solid, blue) 
and one without (dashed, red) the quark mass included. It is clear that the quark mass has 
a large influence on the cross section (the area under the curves). 

In the case of the process gq — > ZQQq, which involves the splitting g —>■ QQ in the 
initial state [Fig. [2](b)] as well as the final state [Fig. [2](a)], this requires that we abandon the 
dipole subtraction method [1] in favor of subtracting the mass singularity via the truncated 
Q distribution function E] 
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Figure 2: Diagrams contributing to the NLO correction to the associated production of a Z 
boson and two high-p^ jets, one or more of which contains a heavy quark {Q = c,b). 

where Pqg{z) = ^[z'^ + (1 — z)"^] is the DGLAP sphtting function. The counterterm is 
constructed by calculating Qq — > ZQq; this cancels the initial-state logarithmic dependence 
on the heavy-quark mass in gq —>■ ZQQq, and yields a cross section in the MS factorization 
scheme. This technique was developed in Ref. 

We have revisited our NLO calculation of Z-l-2 jets, with one or more heavy quarks, and 
corrected this omission. Specifically, we recalculated the NLO processed 

• qq ^ ZQQg at tree level [Fig. [1] 

• gq —>■ ZQQq at tree level [Fig. [2] 

using a finite heavy-quark mass throughout (m^ = 1.4 GeV, mt, = 4.75 GeV). 

The results for the exclusive and inclusive cross sections at the Fermilab Tevatron and 
the CERN Large Hadron Collider (LHC) are presented in Tables II-V. These should replace 
the corresponding tables in Ref. p!]. While there are small changes throughout the tables, 
the only significant change is in the NLO cross section for Zbj at the Tevatron in Table [2] 
(as well as the inclusive cross section Zbj + X in Table H]), which was reduced by 16% by 
including the 6-quark mass (5.23 pb 4.37 pb). The dominant effect of the 6-quark mass 
was on the subprocess qq — ZQQg [Fig. [1], which was reduced by more than 50% (1.16 pb 
— > 0.50 pb). A similar percentage reduction occurs at the LHC (9 pb — > 4 pb), but this 
subprocess is so small at that machine that it hardly affects the total result for Zbj. 

We raised the px threshold on the acceptance for b jets to greater than 15 GeV, and 
found that the effect of the b mass became less. This is consistent with our findings that the 
effect of the c mass is smaller than that of the b mass for the same pt threshold. 

■^It is not necessary to recalculate gg — » ZQQg, since there is no contribution from a virtual gluon splitting 
to a heavy-quark pair, g* — > QQ. We recalculated it nevertheless, and found that the effect of the heavy- 
quark mass is about —6% at the Tevatron and —12% at the LHC. However, since this subprocess makes a 
relatively small contribution to the cross section, this change is negligible. 
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Figure 3: The differential cross section for Z plus two jets, one of which contains a b quark, 
vs. the invariant mass of the b quark and another b quark that lies outside the acceptance 
of the detector (taken to be pt > 15 GeV, \r]\ < 2, ARjj > 0.7), at the Fermilab Tevatron 
(v^ = 1.96 TeV pp). Only the contribution from the subprocess qq Zbbg is shown. The 
(solid, blue) curve includes the b quark mass, while the (dashed, red) curve does not. 
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Table 2: Cross sections (pb) for Z-boson plus two (or more) jets, one or more of which 
contains a heavy quark, at the Tevatron {^/s = 1.96 TeV pp). A jet lies in the range pt > 15 
GeV and \7]\ < 2. Two final-state partons are merged into a single jet if ^Rjj < 0.7. No 
branching ratios or tagging efficiencies are included. Numbers in parenthesis are leading- 
order results. The labels on the columns have the following meaning: ZQj = exactly two 
jets, one of which contains a heavy quark; ZQQ = exactly two jets, both of which contain 
a heavy quark; Z{QQ)j = exactly two jets, one of which contains a heavy-quark pair; 
ZQQj = exactly three jets, two of which contain a heavy quark; ZQjj = exactly three jets, 
one of which contains a heavy quark. For the last set of processes, the labels mean: Zjj = 
exactly two jets, including heavy quarks; Zjjj = exactly three jets, including heavy quarks. 
For ZQj and ZQQ, both the leading-order (in parentheses) and next-to-leading-order cross 
sections are given. The CTEQ6M parton distribution functions are used throughout, except 
for the LO cross sections in parentheses, where CTEQ6L1 is used The factorization and 
renormalization scales are chosen as fip = fiR = Mz- 
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Table 3: Same as Table [2], except at the LHC {^/s = 14 TeV pp). A jet lies in the range 
Pt > 15 GeV and \ri\ < 2.5. (This table is nearly identical to Table III of Ref. [1].) 



LHC 


a (pb) 




ZQj ZQQ ZiQQ)j ZQQj ZQjj 


bottom 
charm 


(352) 418 (109) 92.1 23.5 60.8 92.1 
(443) 624 (87) 75.2 58.6 49.3 123 




Zjj Zjjj 


Z+jets 


(6090) 4840 1810 



Table 4: Inclusive cross sections (pb) for Z boson plus one or two jets, one or two of which 
contains a heavy quark, at the Tevatron (-y/i = 1.96 TeV pp). A jet lies in the range px > 15 
GeV and |?7| < 2. Two final-state partons are merged into a single jet if ^Rjj < 0.7. No 
branching ratios or tagging efficiencies are included. Numbers in parenthesis are leading- 
order results. The labels on the columns have the following meaning: ZQ -|- X = at least 
one jet, at least one of which contains a heavy quark; Z{QQ) = one jet which contains 
a heavy-quark pair; ZQj + X = at least two jets, one of which contains a heavy quark; 
ZQQ + X = aX least two jets, two of which contain a heavy quark; Z{QQ)j = two jets, 
one of which contains a heavy-quark pair. The last row gives the inclusive cross section for 
jets containing both light and heavy partons. The CTEQ6M parton distribution functions 
are used throughout, except for the LO cross sections in parentheses, where CTEQ6L1 is 
used [7]. The factorization and renormalization scales are chosen as = = Mz- The 
uncertainties are from the variation of the renormalization scale, the factorization scale, and 
the parton distribution functions, respectively. 



Tevatron 


a (pb) 




Z + 1 jet +X 


Z + 2 jets +X 


ZQ + X Z{QQ) 


ZQ] + X ZQQ + X Z{QQ)] 


bottom 
charm 


(8.23) 18.1 2.1 
(11.3) 27.5 6.6 


(0 A 7n+l-2 +0.07 +0.5 (0 A7\ "i. 7/1+0.45 +0.12 +0.15 n fiQ+0-26 +0.06 +0.05 
I^Z.iOJ ' U_o.9 -0.05 -0.5 l^-^'J -J- '^-0.45 -0.12 -0.15 U.O-3_o.l6 -0.06 -0.05 

(''^ On^l 7 8e;+l-7+0.15 +0.8 /o o/i \ q q7+0.44 +0.0 +0.14 9 nn+0-85 +0.22 +0.15 
\O.Z,\J) ( •O0_;^ 4 _Q 26 -0.8 l^-'J^J 'J-'J ' -0.37 -0.87 -0.14 ^-UU.q.ss _0.16 -0.15 


all jets 


(898) 1070 


(163) 2Qht'^tltl 
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Table 5: Same as Table HI except at the LHC (a/s = 14 TeV pp). A jet lies in the range 
Pt > 15 GeV and |?7| < 2.5. (This table is nearly identical to Table V of Ref. [1].) 



LHC 


cr(pb) 




Z + 1 jet +X 


Z + 2 jets +X 


ZQ + X Z{QQ) 


ZQj + X ZQQ + X Z{QQ)] 


bottom 
charm 


(826) 1060 25 
(989) 1430 50 


(352) 510i^t ni (109) 1531^° tl 24t^° t\i 
(443) 747tl}° t% til (87) 125t\l tl tl 591^ tl tl 


all jets 


(15300) 18400 


(6090) QQ^Otll'^oni'ntl 
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